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ABSTRACT 

Do molecular clouds collapse to form stars at the same rate in all environments? In large spiral galaxies, the rate of transformation of 
H2 into stars (hereafter SFE) varies little. However, the SFE in distant objects (z ~ 1) is much higher than in the large spiral disks that 
dominate the local universe. Some small local group galaxies share at least some of the characteristics of intermediate-redshift objects, 
such as size or color Recent work has suggested that the Star Formation Efficiency (SFE, defined as the SFRate per unit H2) in local 
Dwarf galaxies may be as high as in the distant objects. A fundamental difficulty in these studies is the independent measure of the H2 
mass in metal-deficient environments. At 490 kpc, NGC 6822 is an excellent choice for this study; it has been mapped in the C0(2-l) 
line using the multibeam receiver HERA on the 30 meter IRAM telescope, yielding the largest sample of giant molecular clouds 
(GMCs) in this galaxy. Despite the much lower metallicity, we find no clear difference in the properties of the GMCs in NGC 6822 
and those in the Milky Way except lower CO luminosities for a given mass. Several independent methods indicate that the total H2 
mass in NGC 6822 is about 5 x 10'' Mq in the area we mapped and less than 10^ Mq in the whole galaxy. This corresponds to a 
A'(H2)//co ~ 4 X 10^' cm"^/(Kkm s"') over large scales, such as would be observed in distant objects, and half that in individual 
GMCs. No evidence was found for Hi without CO emission. Our simulations of the radiative transfer in clouds are entirely compatible 
with these A'(H2)//co values. The SFE implied is a factor 5-10 higher than what is observed in large local universe spirals. The CO 
observations presented here also provide a high-resolution datacube (1500 a.u. for the assumed 100 pc distance, 0.41 km s"' velocity 
resolution) of a local molecular cloud along the line of sight. 
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1. Introduction 
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In the study of star formation in a cosmological or extragalactic 
context, rather than the details of the collapse of a cloud core to 
a star, we are interested in understanding why stars form where 
, they do, whether the efficiency varies, and what factors influ- 
' ence the initial mass function (IMF) of the stars. Over the last 
I decade, it has become very clear that the star formation rate 
per co-moving volume was much higher in the past, some 10 
or 20 times the current rate at a redshift of z ~ 0-5..~ ^ ( ^-S- 
; lMadauetal.lll996t iHeavens et aLll2004t IWilkins et al.ll2008h . In 
I turn, this shows that the transformation rate of gas into stars was 
considerably (factor few at least) higher when the universe was 
roughly half its current age. Galaxies at that time were smaller 
and of lower metallicity, such that naively at least one would ex- 
pect that th e molecular-to-atomic gas mass ratio would be lower 
than today dYoung & Knezeklll989HCasoli et al.lll998h . making 
the higher efficiency even more surprising. Since stars form from 
H2, and not directly from H I (with the possible exception of the 
so-called Pop. Ill, or first generation stars), this suggests that ei- 
ther large amounts of molecular hydrogen were available or that 
for some reason the efficiency of star formation (SFE, defined 
as the SFR per unit H2) was particularly high back then. In fact, 
the SFRs proposed are so much higher than the SFR today that 
both possibilities may be required. Because at least 10% of the 



baryons in galaxies today are thought to be in neutral gas (and 
more than 10% in many cases), an SFR a factor 15-20 higher 
must result at least partially from a higher SFE. If the SFE is 
higher, then something about the process of star formation is 
different and there could be other important differences like a 
change in IMF. Moderate to high redshift galaxies are typically 
smaller and more gas-rich than today's spirals and most likely 
have a slightly subsolar metallicity. They thus resemble today's 
small spirals such as M33, or the even smaller NGC 6822, and 
could be expected to have a low H2/HI mass ratio. If so, this 
would make the SFE in these objects even more extreme. 

The first step is to learn more about the molecular gas content 
of galaxies with these properties. Significant quantities of molec- 
ular gas were d etected far out in the outer disks of NGC 4414 
and NGC 6946 dBraine & Herpinll200l iBraine et aI.lF2007i) . The 
outskirts of spirals share the subsolar metallicities and low mass 
surface densities of small and/or medium/high redshift spirals 
but not the level of star formation. We are fortunate to have a 
number of small galaxies in the Local Group close enough that 
individual giant molecular clouds (GMCs) can be resolved and 
without distance ambiguities. A first step has been taken, show- 
ing that molecular gas forms very far out in M 33 despite the low 
metallicity and very low ambient pressure (G_ardan et al. 200J^. 
NGC 6822 is among the nearest galaxies and i s a small late -type 
dwarf spiral at a distance of about 490 kpc (lMated[T998h and 
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has a mass and luminosity of roughly 1 % of that of our galaxy, 
thus representing a step down in mass, luminosity and metal- 
Ucity (roughly 12 + log(0/H) ~ 8.1 (Lee et al. 2006)) with re- 
spect to M33, itself a step down in the same quantities from the 
Milky Way or M31. At 490 kpc, 1" corresponds to 2.4 pc, such 
that Giant Molecular Clouds (GMCs) can be resolved with large 
single-dish radiotelescopes. 

Observations of the universe at redshifts z ^ 1 show that to- 
days large spirals were not present or rare at these earlier epochs. 
Rather, the galaxies were smaller and had higher star forma- 
tion rates. They probably had somewhat lower me tallicities and 
were bluer. M33 was found to have a high SFE (iGardan et alj 
l2007h . compared to local un iverse spirals and ICIO appears to 
show a high SFE as well dLerov et al.l l2006h . NGC 6822 is 
closer and its molecular gas content has not been mapped sys- 
tematically until now. Giant Molecular Clouds have been re- 
solved and their physical properties studied in Local Group 
galax ies both w ith single dish tele scopes (M 33 by Gardan et al.l 
120071: SMC by iRubio et al.lll993t LM C by iFukui et al.ll2008F 
SMC and LMC by llsrael et al."2003b|) and with i nterferome- 
ters (M3 3 by .Engargiola et al. 2003; IC 10 bv lLerov et al. 200^ 
M31 bv iRosolowskvl 120071) . NGC 6822 has been observed at 
many wavelengths to study the interstellar medium (ISM), the 
dynamics, and trace the star f ormation. Spitzer FIR observations 
were carried out recently by Cannon et al. ('2006'). The atomic 
gas has been mappe d (de Blok & Walter 2000, 2003, 2006a^ 
IWeldrake et ai]|2003l) and the molecular gas observe d at specific 
positions with the 15m SEST and JCM T telescope (llsrael et al.l 
[l996t [Israeli 1997at llsrael et al.ll2003ah .With the OVRO interfer- 
ometer (IWilsonlll994 observed 3 GMCs in NGC 6822 in the 
Hub V region. Since one of the main questions is whether sig- 
nificant quantities of molecular gas could be present without de- 
tectable CO emission, high sensitivity high resolution mapping 
of large regions allowing the detections of individual possibly 
optically thick clouds is required. The molecular gas content de- 
rived via CO can then be compared with other means of tracing 
the molecular and atomic gas. 

In this article we present the observations and data reduc- 
tion, mostly of CO (but also '^CO and HCN(l-O) in Hubble V), 
followed by the production of a catalog of molecular clouds and 
their properties which we compare with Galactic GMCs. Two 
methods were used to compile the catalog of cloud sizes, CO 
intensities, and virial mass es: vis ual inspection of the data cube 
and the CPROPS algorithm ( iRoso lowskv & Lerov 2006). A map 
of the total CO emission is then compared (Sect. |6]l with other 
means of estimating the H2 column density, leading to a discus- 
sion of the SFE in N GC 6822. We then p resent several models 
usmg the CLOUDY (iFerland et al.ll 19981) code to compare sev- 
eral spectral synthesis models with the observations presented 
in the preceding sections. Finally, two regions of NGC 6822 are 
discussed in more detail - Hubble V and X. A local molecular 
cloud is present along the line of sight to NGC 6822 and has 
thus been observed serendipitously at high spectral and spatial 
resolution. These data are presented in the last section. 

2. Observations 

NGC 6822 was observed during three separate runs at the IRAM 
30 meter telescope in November 2006, February/March 2007 
and August 2008 in mostly good weather All mapping was 
done using the HERA array of 9 dual polarization receivers 
(Schuster et al. 2004) in the CO(2-l) line, whose rest frequency 
is 230.53799 GHz which gives a nominal resolution of 11.4" 
for this line. The On-The-Fly mode was used to cover a roughly 



Table 1. Parameters for NGC 6822. 



Q'<,(J2000)" 
(5„ (12000)" 
Distance'' 
Vlsr 

Metallicity'' 

Total optical luminosity'' 
Total H I mass'' 
Dynamical mass*^ 



19i'44"'57 ^:83 
-14°48'11" 
490 ± 40 kpc 
-47 kms"' 
^ 0.3 Zq 
~ 9.4 X 10'' Lq 
1.34 X 10** Mq 
3.2 X 10" Mq 



(0,0) o ifset of our HERA CO(2-l) map 

Mateo (T99I) 

Skillman et"al| (il989l): lLee eTaTI j2006l) 

de Blok & Walter (2066b) 

de Blok & Walter, ( .2006h) : IWeldrake et alj ( 120031) 



11x7 arcminute region, scanning along the RA and then Dec di- 
rections and observing a reference position offset from our cen- 
tral position by (100", 500") before and after each scan. The ref- 
erence position was chosen to be outside of the 5 x 10^" cm"^ 
contour of the H I column density map (see Fig.[T|i. 

The VESPA backend was used with a channel spac- 
ing of 312 kHz or 0.406 km s ' covering velocities from 
-1-20 to -150 km s ', well beyond the rotation curve of 
NGC 6822. Local (Galactic) emission was detected around 
-1-5 km s"' (see Sect.lTOli. All data are presented in the main beam 
temperature scale and we have assumed forward and main beam 
efficiencies of rjfor - 0.90 and ^,„;, = .52 for the HERA ob ser- 
vations, the sensitivity is then 9.6Jy/K (ISchuster et al.ll2004l) . 

Also during the November 2006 run, during poorer weather 
than for the more demanding HERA observations, the major 
HII region Hubble V was observed in the '^00(1-0), '^^CO(2- 
1), HCN(l-O), and '^00(1-0) fines. Wobbler switching was 
used with a throw of 60" and the 100 kHz and VESPA back- 
ends were used, yielding spectral resolutions of respectively 
0.27, 0.43, 1.06, and 0.26 km s"' for the fines above. Tfie for- 
ward and main beam efficiencies at these frequencies are as- 
sumed to be respectively rjfor = 0.95, 0.91, 0.95, 0.95 and 
T]„,h = 0.75, 0.55, 0.78, 0.74. The data reduction of the HERA 
observations is described in the next section. For the Hubble V 
data, bad channels were eliminated and spectra were averaged, 
yielding tfie spectra discused in Sect. 19.2! 



3. Reduction of HERA data 

Tfie On-Tfie-Fly mapping tecfinique witfi a multi-beam array 
generates a fiuge amount of spectra, more tfian a million in tfie 
case of NGC 6822. Inspection of individual spectra is tfius not 
possible and tfie reduction was automated. All data reduction 
was done witfiin tfie Gilda^ CLASS and GREG software pack- 
ages. After filtering out tfie spectra taken in very poor condi- 
tions (Tjvs > 1500 K), we treated tfie main problem wfiicfi was 
tfie sligfit platforming wfiere tfie sub-bands of tfie auto-correlator 
backend were stitcfied togetfier. Tfie platforming effect is tfie re- 
sult of different non-linearities in tfie sampling stages of tfie sub- 
bands. In case of cfianging total power levels as compared to tfie 
reference position tfiis introduces offsets in tfie subbands. Tfie 
steps were very small in our case but sufficient to affect weak 
lines. For eacfi spectrum, tfie average value for eacfi subband, 
outside of tfie line windows as far as possible, was subtracted 
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-14°40'00" 



-15°00'00" 



Fig. 1. Main image GALEX FUV filter image of NGC 6822. The yellow box indicates the are a observed in CO. The das hed, dotted, 
and solid yellow contours indicate the HI column density at 12" x 42" resolution fromlde Blok & Waited (l2006bl) at column 
densities of 0.5, 1, and 1.5 x 10^' cm"^ respectively. The white cross corresponds to (0,0) offsets of the HERA CO(2-l) map. 
Bottom right HI contour at a column density of 0.5 x 10~' cm"^. The dashed rectangle corresponds to the GALEX FUV map above, 
the solid line contour to the area observed in CO. 



from each channel of the sub-band, eliminating the platforming. 
This process takes out a zero-order baseline. 

There is a known baseline ripple with HERA dSchuster et alj 

|2004|) corresponding to a reflection off the secondary mirror at 
6.9 MHz or 9 km s Since the pixels (receivers) are affected at 
quite different levels and because the width of the ripple is close 
to that of molecular clouds, we tested each spectrum and when 
the ripple was strong enough to be identified in the individual 
(i.e. roughly Isec integration time) spectra, the corresponding 
Fourier frequency was replaced by an interpolation based on ad- 
jacent frequencies. This fairly standard filtering is often applied 
"blindly" to all spectra but we only applied it when necessary. 

In order to create the datacube, we created a data table 
with the TABLE command within CLASS90 and then used the 
XY JVIAP task in GREG with parameters such that the final res- 
olution became 15", or about 36 pc at the distance of NGC 6822. 
Data cubes with different resolutions where always generated di- 



rectly from the original data by convolution with the correspond- 
ing kernel size. 

4. The individual molecular clouds in NGC 6822 

4.1. Identification by visual inspection 

The final beam width of the CO observations is about 30 pc and 
because of the beam dilution we do not expect to see clouds 
with sizes under 10 pc. We define clouds as structures similar to 
Galactic Giant Molecular Clouds (GMCs) that appear as grav- 
itationally bound and non transient structures a few ten of par- 
sees in size. Figure |2] shows the CO(2-l) and HI spectra of the 
clouds found in the original 15" datacube, showing both the CO 
and HI intensity scales. The spectra are averages over the 50% 
brightness contour It is clear that all strong CO lines are close 
to the HI peak in velocity. However, cloud 13 is at either the 
edge of the H I line or possibly part of a second H I feature with 
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-14°47'30" - 



-14°47'40" - 



-14°47'50" - 




-14°48'00" - 



Fig. 3. Contour integrated intensity maps of clouds 4 (solid lines) 
and 5 (dashed lines) contours every O.IK km s"'. The emission 
peaks are separated by about 9 arcsec showing that there are 
clearly two separate clouds. 



a brightness of 25 K. Cloud 14 is also near the edge of the HI 
line although mostly within the 30 K H I brightness temperature 
level used to define the CO line window in the next section. 

Gaussian fits were made to the individual clouds in order 
to determine linewidths and central velocities. In several cases, 
more than one gaussian was required. Since one of the goals 
is to measure line width and cloud sizes for NGC 6822 clouds 
and compare to Milky Way clouds, only the stronger (in antenna 
temperature) and narrower gaussian was used to define the line 
widths and sizes of the individual clouds. Line areas are com- 
puted by summing channels in a velocity range determined man- 
ually for each cloud. The line intensities, systemic velocities and 
widths are averages over the 50% brightness level of each cloud. 

Spectra for clouds 4 and 5 show that these clouds are only 
partially spatially resolved. Contour maps of the integrated in- 
tensity for the two components 4 and 5 (Fig.|3]l indicate that the 
emission can be separated into two clouds separated by about 9 
arcseconds. 

Table 2 provides positions, velocities, linewidths, and esti- 
mated sizes for each cloud. In some cases, clouds were sep- 
arated by summing over different velocity ranges, enabling 
cloud separations smaller than the beamsize. On the individ- 
ual maps, corresponding to each cloud, the 50% brightness 
level was defined for each cloud, allowing the effective ra- 
dius (see Sect. 14.2.1b of the cloud to be estimated as Re - 

(1.92/2 V2 ln2) ^ polygon area" - 15"^ and converted to pc. 

The subtraction of the beam area (15~) enables a simple decon- 
volution with the the beamsize assuming in average a gaussian 
intensity distribution of the cloud emission, the numerical fac- 
tor converts the FWHM to an effective radius. Below about 10 
pc, the beam dilution makes the detection of individual clouds 
difficult. 

The CO(2-l) integrated intensities (col. 7) are obtained by 
summing channels within V - Vg ± AV. Virial masses (col. 

8) are calculated from M „,v = 1040/?, (Ay/2 V21n2)' follow- 
ing ISolomon et al.l (Il987h and using their form factor of 2.7. 



For comparison, cloud masses (including He) can be estimated 
as M = /co(2-i) X 2500 x (size/2Opc)^M0 using a "standard" 
Milky Way factor of jy (H9)//mri-n^ = 2x 102°cm-2/(K km s"') 
(iDickman etaLl 119861) and a CO(2-1)/(1-0) Une ratio of 0.7 
(ISawada et al.l 120011) . Inspection of the values shows that the 
Galactic A^(H2)//co factor yields masses far below the virial 
mass es. This is well known for clouds in low metallicity galaxies 
re.g. iRubio et al.lll993l: Ilsraellll997bh . The virial masses them- 
selves may be underestimates of the true H2 masses if the H2 
extends beyond where the CO is detected. The following col- 
umn provides an independent estimate of the H2 mass (see Sect. 
l6.3l for details) and resulting N(H2)/Ico factor. 



4.2. Automated cloud identification 



4.2.1. CPROPS 



We have also used the CPROP^ program dRosolowskv & Lerovl 
2006) to identify GMCs and measure their physical properties 
in an unbiased way. CPROPS first assigns the measured emis- 
sion to clouds by identifying emission above a 4cr noise level 
and then decomposing the e mission into individual clouds (see 
iRosolowskv & Lerovll2006l for further details). It then extrap- 
olates physical properties such as cloud sizes and masses to 
OK noise level, independent of the beamsize (i.e. deconvolved). 
Using CPROPS, we find 1 1 clouds all of which have also been 
identified as such by eye. The cloud properties as identified by 
CPROPS are presented in the lower part of Table |2] We were not 
able to setup CPROPS in such a way that all the eye-identified 
clouds were found without a large number of clouds we do not 
believe are real being also identified by CPROPS. CPROPS was 
used with the following parameters: a constant distance DIST 
equal to 490 kpc; the /NONUNIFORM parameter along with a cus- 
tom noise map computed from velocity channels without any 
signal (i.e. outside the rotation curve), to take into account non 
uniform noise over the CO map; the following values for the de- 
composition parameters FSCALE=2.0, SIGDISCONT=Q to en- 
sure that each area of non-contiguous emission is assigned to an 
individual independent cloud. 

Using the identified cloud emissions, CPROPS computes four 
initial quantities through an extrapolation to a K noise level: 
the cr,- and cr,, spatial dispersions, the velocity dispersion cr,, and 
the CO flux Fco- From these quantities, the following physical 
quantities are deduced. An effective radius (col. 6) is obtained 
using the formula: 



(1) 



iFWHMbear. 

\ 2V21n2 



i FWHMhear 

\ 2V21n2 



(2) 



Where 77 - 3 .4 / ^J^i - 1 .92 is an empirical geometric factor from 
ISolomon et al.l (Il987h to take into account the radial distribution 
of the gas density inside the molecular cloud and cr„„„ and cr„„j 
are spatial dispersions along the major and minor axes of the 
cloud. 

For a few of the clouds identified by CPROPS, the minimum 
dispersion cr„„„ was found to be smaller than the beamsize. In 
these cases we have chosen to use an arbitrary value of cr,„,>, cor- 
responding to an effective minimum spatial dispersion of lOpc. 



|http : //www ■ cf a . harvard . edu/- erosolow/CPROPS/| 
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We have also computed an extrapolated radius (col. 7) 

jA2a - (FWHMbea,n/2 V2lK2)2 

Re-c = (3) 

7T 

where A20- is the area of the individual clouds extrapo- 
lated down to a 2cr level projected on the sky plane (see 
iRosolowskv & Lerovll20 06. for details). 

The clouds' integrated CO(2-l) intensities (col. 8, lower part 
of Table|2]i were obtained for each cloud by dividing the CPROPS 
computed Leo luminosity by the projected area of each individ- 
ual cloud, A21T, at the two sigma level. 

4.2.2. Comparing with eye identification 

For the subset of clouds which have both been identified by eye 
and by the CPROPS package, we can compare the physical prop- 
erties obtained by the two independent methods. Table |2] shows 
the properties computed by hand for the 15 clouds identified by 
eye, and by CPROPS for the first eleven which have also been 
identified with CPROPS. Each property is computed slightly dif- 
ferently for the two methods and the next paragraphs will explain 
these differences. 

In the manual identification method, the offsets for the cloud 
positions are obtained by taking the average position in right as- 
cension and declination of the pixels inside of the half maximum 
contour of the cloud emission. In the case of CPROPS, the offsets 
are equal to the first moments of the cloud emission (down to 2cr) 
along the right ascension and declination axis. The positions are 
in general the same to within 0.2 beam fwhm. 

The systemic velocity is taken as the average of the gaus- 
sian fit to the CO line (see Sect. 14.1b for the manual identifica- 
tion method, in the case of CPROPS it is computed as the first 



moment along the velocity axis of the cloud emission down to 
2cr. Following the same idea, the line width are computed in the 
manual case as the the function width at half maximum of the fit- 
ted gaussian and as the second moment converted to FWHM by 
multiplying by 2 V21n2 for the properties derived from CPROPS. 
An effective radius Re is obtained using Eq.|2]and preceding sec- 
tion. In the case of the manual identification, the effective radius 
was obtained in the following way. 

1 92 ; 

R?' = ; J A I - Ahean, (4) 

2V21n2 > - 

where A ^ is the area in pc^ inside the contour at half of the peak 

integrated intensity for each cloud. And Abeam - ITlOpc^ is the 
area subtended by the 15" beam at a distance of 49Qkpc. The 
2 V2 In 2 factor converts the FWHM to dispersion and the 1 .92 
from dispersion to effective radius. 

The CO intensity was computed using CPROPS CO luminos- 
ity and dividing it by the area A20- of the cloud computed above. 
In the case of the manual selection, the emission inside the 50% 
level contour was summed over the velocity range V,, + AV and 
multiplied by the 0.41 km s ' channel width. 

In both cases, the virial mass was ob tained using the follow- 
ing formula from dSolomon et alJ[T987h : 

/ Ay \^ 

M,,,.^ 10407?, ^ (5) 

\2 V21n2/ 

The uncertainties in the virial mass estimates are dominated by 
the hypothesis that the molecular clouds are indeed gravitation- 
ally bound and by the value of the geometric factor describing 
the density distribution of the gas. The marginally gravitation- 
ally bound case of a cloud in isolation with no magnetic field 
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Table 2. Properties for the clouds identified by eye (Top), and by the CPROPS package (Bottom). 
The clouds whose numbers are in bold face correspond to our final sample of clouds 



Cloud 




<5off 


Vlsr 


FWHM 


Re 


Ico 


Mvi,-" 


Mh/- 


A'(H2)/Ico'' 




arsec 


arcsec 


km s"' 


km s"' 


pc 


Kkms-' 


Mo 


Mo 


K km s-1 



GMC NGC6822 1 


-48 


330 


-43.2 


1 


.5h 


tO, 


.9 


22!?4 


0, 


Q9+0.10 


Q ^+24 


X 


10^ 


3 6+'* 


X 


10" 


1 4+0.2 
-^•^-0.2 


X 


10^' 


GMC NGC6822 2= 


-68.5 


299 


-42.4 


5, 


.1 H 


tO, 


.4 


291 


2, 


47+0.14 

•^'-0.13 


1 9+"^ 


X 


10' 


3 1+"' 


X 


10" 


3 0+" ' 


X 


10^° 


GMC NGC6822 3 


80 


58 


-38.0 


2, 


.4h 


tO, 


.5 


38!^^ 


0, 


35+0.07 


4.1+f 


X 


10-* 


'■-'-3.7 


X 


10" 


4 3+"= 
^■-'-1.1 


X 


10^' 


GMC NGC68224 


39 


25 


-39.9 


2, 


.Oh 


tO, 


.6 




1 


34+0.10 


1 2+"' 

i.z. „3 


X 


10' 


1 2+"-' 


X 


10' 


■2 1 +0.4 
-'■^-0.3 


X 


10^1 


GMC NGC6822 5 


35 


24 


-43.2 


3, 


.4h 


t 1, 


.1 




1 


<;9+0.11 


4 1 +0.6 
^■^^-0.3 


X 


lO'' 


1 2+''-2 


X 


10' 


^■-'-0.2 


X 


10^' 


GMC NGC6822 6 


30 


16 


-48.3 


1 


.6h 


tO, 


.4 


12t'.'. 


0, 


09+0.05 
•"^-0.07 


5.8+.'.''. 


X 


10^ 


4 7+' "' 
^■'-1.4 


X 


10" 


9 4+0.3 


X 


102' 


GMC NGC6822 7 


-10 


364 


-38.1 


3, 


.2h 


tO, 


.4 


27+11 


0, 


zLS+0.06 


5 2+-* ' 
-'•^-34 


X 


lo-* 


4 Q+0.2 


X 


10^ 


2.8+|-^' 


X 


102« 


GMC NGC6822 8 


31 


61 


-44.8 


2, 


.4h 


1 1 


.4 


14+28^ 





CQ+0.13 
•-"'-0.14 


1.5+.'.''. 


X 


10-* 


4 1+5' 


X 


10" 


2-8!i;» 


X 


102' 
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" Offsets with respect t o the reference positio n (a„, 60) = (19''44'"57 ';8, -14°48' 1 1") 
'' Virial mass following lSolomon et al.l l fl987h using Re and AV (see Sect.|4ll 

H2 mass estimate from 8/jm emission (see Sect. 16.3b . 
'' X conversion factor for individual clouds using A'(H2) estimated from 8/jm emission (see Sect. l6.3l >. 
= Hubble V( Hubble 192^) 
f see Sect.lH 

8 The CPROPS computed FWHM is clearly overestimated compared to the one computed by hand by a factor of about 3. Myjp being proportional 

to FWHM", this explains the order of magnitude difference between the computed virial masses for this cloud. 
'' The computed H2 mass for these clouds was found to be negative. 
' The deconvolved size was not defined. 



would yield masses a factor 2 less than virial. The virial masses 
are widely used because clouds have magnetic fields, are not iso- 
lated, and collapse to form stars. 



The H2 masses (col. 9 for top and col. 10 for bottom parts of 
Table|2]i are derived following the method described in Sect. 16.31 
using the 50% contour in the case of the manual identification 
and the 2cr contour for CPROPS. 



The last column in Tab[2l A^(H2)//co is the H2 column den- 
sity derived in Sect. l6.3l divided by Icq. 



4.3. Estimation of errors 

For the CPROPS identification method, the uncertainties for 
the quantities FWHM, Re, Ico and Mvii are computed us- 
ing the bootstrapping method of CPROPS (see details in 
Rosolowsky & Lerov 2006). We now describe the computation 
of the uncertainties in the case of the identification of clouds by 
eye. The FWHM line width and associated error are computed 
with the CLASS software gaussian line fitting algorithm. Then, 
using the line widths obtained, the unvertainty in integrated in- 
tensity cTKkmis is calculated for each cloud. The uncertainty on 
the size Rg is estimated by calculating Re from contours (cloud 
sizes) placed at Ip<,ak/2-o-Kkm/., and Ipeakl^ + CTKkmis, thus brack- 
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Fig. 4. Size vs linewidth for molecular clouds in NGC 6822, 
open and filled circles correspond respectively to measured 
and deconvolved sizes (see Sect. |4|. Filled squares correspond 
to physical properties computed with the CPROPS package. 
AV PC Rg is the galactic relationship from ISolomon et al.l 
(Il987h . Notwithstanding considerable scatter, the distribution of 
the clouds (corrected for finite beamsize) in this diagram appears 
consistent with a size-Unewidth relation similar to that in the 
Milky Way. 




Lco(2-i) (K km/s pc ) 

Fig. 5. Virial mass vs CO luminosity for molecular clouds in 
NGC 6822, open and filled circles correspond respectively to 
measured and deconvolved sizes (see Sect. |4]i. Filled squares 
correspond to physical properties computed with the CPROPS 

package. The Une Myj^ ^ 9796 ( ^^tmt ) 

is the galactic relation- 
ship from lSolomon et aTl (Il987h . the dashed line corresponds to 
a constant galactic A^(H2)//co ratio of 2 x 10^" cm"^/(K km s"'). 
The virial masses are systematically larger for a given luminos- 
ity in NGC 6822 compared to clouds in the Milky Way. 



eting the cloud size obtained using the I peak I ^ contour. This gives 
respectively an upper (lower) bound on the value of Rg. The er- 
rors are then propagated into the Virial mass. For each cloud, 
the contours defined at I peak 1'^ ± o-Kkm/s are used to compute the 
molecular gas mass (see Sect. 6) from the 8 micron map and thus 
estimate the uncertainties in the molecular gas mass (Tab. 2 col. 
9) and the A^(H2)//co factor (Tab. 2 col. 10). 

5. The size-linewidth relation for the molecular 
clouds in NGC 6822 

Figures |4] and |5] show respectively the size-linewidth relation 
(Ay vs. Re) and the virial mass vs CO l uminosity, showing in 
both cases the Galactic values taken from lSolomon et all (Il987l) 
as a straight line. The distribution of the clouds in NGC 6822 
appears consistent with a size-linewidth relation similar to that 
in the Milky Way GMCs. In recent work lHever et all (l2009h ob- 
tain lower H2 masses and a dependency on the square root of 
the surface density, the variation we obta in in these parameter i s 
much to low to reach a conclusion. Like lSolomon et al] ( Il987h . 
they concluded that GMCs are gravitationally bound. 

Fig.|5]shows that there is a factor several difference between 
the virial masses of the NGC 6822 clouds and the masses ob- 
tained from the size and CO(2-l) intensities using a Galactic 
A^(H2)//co conversion factor For a CO(2 - l)/CO(l - 0) ratio 
of 0.7, the diff'erence is 5 - 6, the true N{H2)/Ico value is then 
at least 5-6 times higher than in the Milky Way. Due to the low 
metallicity of NGC 6822, cloud sizes as seen in CO are probably 
underestimated compared to Galactic observations, where the 



shielding will be much more efficient to protect CO molecules 
and allow the CO size to be similar to the total size of the H2 
dominated region (the molecular cloud). Because the line width 
is presumably determined by the total mass, the virial masses 
should be underestimated linearly with the size. CO luminosi- 
ties, used when applying a N{H2)/Ico factor, should be underes- 
timated twice as much in proportion because luminosities vary 
with the square of the size. In the following section we try to 
estimate the A^(H2)//co(i-0) factor by means without these draw- 
backs. 

6. The total molecular mass of NGC 6822 

In this section we describe how we make a CO integrated inten- 
sity map to trace the H2 column density. In order to test whether 
there could be substantial amounts of molecular gas without 
associated CO emission , we use two other alternative methods 
(similar to lIsraeliri997al) to estimate H2 masses for comparison. 

6. 1 . CO(2- 1 ) intensity maps 

We compute the CO(2-l) integrated intensity map using a mask- 
ing method, taking into account the 21cm atomic hydrogen line 
data, we developed in order to filter out some of the noise 
present in the observations and increase the sensitivity to low 
intensity possibly diffuse CO emission. Previous masking meth- 
ods used masks created from spatially smoothed versions of the 
original CO data cubes to filter out regions dominated b y noise 
(I Adler et al Jl 1991 iDi gel et aPl 1 996t iLokiard et al.ll 1 999h . 
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Fig. 6. CO(2-l) intensity map obtained using the masking method described in Sect. 16.11 The map unit is K km s , the yellow 
HI contour is at A^(HI) = 10"' cm"^, the red contours are the half maximum levels of each cloud identified by eye, the numbers 
correspond to the ones in Table |2] The cross corresponds to (0,0) offsets of the HERA CO(2-l) map. 



We use the 21cm atomic hydrogen data a t 12" x 42" x 
1.6 km s ' resolution dde Blok & Waited l2006ah to achieve the 
same goal, the underlying hypothesis being that molecular gas 
is unlikely to be present for low enough values of A^(H I) so the 
corresponding velocity channels can be discarded when com- 
puting the integrated intensity CO map. For each pixel of the 
H I cube, we estimate a noise level from velocity channels that 
clearly contain no signal from NGC 6822, we then create a bi- 
nary mask keeping only the velocity range for each pixel corre- 
sponding to a H I signal value above a defined factor of the pixel 
noise. Since the noise in the H I cube varies little over the region 
observed in CO, a cut in S/N is like a cut in antenna tempera- 
ture. The integrated moment map for the CO(2-l) data (Fig.|6]l 
is then computed summing only velocity channels included in 
the H I mask. The result is an increased S/N ratio as the channels 
contributing only noise to the sum are no longer taken into ac- 
count. The value of the noise threshold was chosen at 6cr which 
corresponds a map averaged H I brightness temperature of 30K. 
We tested masking values between 25 and 40 K (5 to 8 cr) and 
the total CO intensity varied by only a few percent. Significantly 
above or below these values, CO signal is lost or more noise is 
included. Using this procedure we miss the very weak cloud 13 
and part of cloud 14 shown in Fig.|2] 

The values in the CO integrated intensity map (Fig.|6]) yield 
H2 column densities when multiplied by a A^(H2)//co factor. If 
we sum all of the emission in Fig.|6l we obtain about Lco(2-i) ~ 
5 .5 X 1 0"* K km s" ' pc^, or some 3x10^ M© in the region we have 



observed for a NiHoJ/Ico ratio of 2.0 x 10^° cm-V(K km s"'). 
Our HERA map covers an area corresponding to 40% of the to- 
tal Spitzer lOfim luminosity of NGC 6822 but over 60% at the 
IMjy/sr cutoff we apply later to be less affected by the noise. 
Thus for all of the galaxy we can expect the total CO(2-l) lu- 
minosity to be between 1.5 and 2.5 times our value. Assuming 
a ratio of 0.7 between CO(2-l) and CO(l-O) we can estimate 
a total CO(l-O) luminosity of 16 ± 4 x 10"^ km s"' pc" in sub- 
stantial agree ment with the value 12(+12, -6) x 10^ km s ' pc^ 
[Israeli (Il997ah estimates for the whole galaxy. The CO emission 
is thus rather weak and the next step is to compare with the other 
means of locating molecular gas. 

6.2. Infrared data 

All of the infrared maps are take n from the SINGS (S pitzer 
Infrared Nearby Galaxies Survey, iKennicutt et al.l l2003h fifth 
public data delivery. The 8jum map does not show a morpho- 
logical similarity with the local galactic emission as seen in CO; 
we have therefore neglected the Milky Way cirrus contribution 
in this band. For the 160yu m , IQ/um and 24^m MIPS band we 
have used the maps from dCannon et al ] l2006h where a smooth 
component fitted on emission outside NGC 6822 representing 
the local emission and instrumental and observational bias has 
been substracted. The IRAS lOOjum data was obtained using the 
IPAC HIRES algorithm using default parameters with no addi- 
tional smooth component subtracted. 



p. Grader et al.: The Molecular ISM of NGC 6822 



9 



N(H2) from Sg^ (cm"') 




19^45'"00' 40' 



Fig. 7. H2 column density map estimated from PAH 8jum emis- 
sion, CO(2-l) intensity contours every (O.lSKkms"') starting 
from (0. 15K km s"' ). The yellow rectangle is the region we have 
mapped in CO(2-l) with HERA. The beams for the NH2, CO are 
plotted in white. 



6.3. Molecular gas mass from Sfim emission 

Before attempting to use the 8 /im emission to trace gas, the 
8/im data was corrected for stellar continuum emission by sub- 
tracting the 3.6iUm emission scaled by a factor 0.232, following 
iHelou etal] ( |2004|) . In the 8//m PAH band, we expect that less 
than 10% of the emission is from hot dust. This comes from ex- 
trapo lating a blackbody c urve from the 24/im point in Fig. 12 
from iDraine '&n (l2007l assuming the dust emitting at 24yum 
can be considered "hot") to 8yum and comparing that to the 8/jm 
emission on the curve. At most, this would slightly reduce the 
gas masses we calculate for Hub V and X. We have therefore 
not subtracted a hot dust contribution from the 8jum PAH emis- 
sion. 

PAHs have often been used as tracers of Star Formation, in- 
cludi ng in high re dshift objects for which very little is known 
(e.g. lAussel et al.|[l999) . In this section we propose to use the 
PAH emission observed in band IRAC4 (Spitzer) to trace the 
gas. There is a close theoretical relationship between the PAH 
band emission per H-atom (atomic or molecular) and the inci- 
dent U V field up to UV fields of 10000 times solar dPraine & Lil 
I2OO7I Fig. 13, lower panel) The UV fields in NGC 6822 are far 
below this value. In cloud cores, where UV emission is not avail- 
able to excite the PAHs, very little 8//m emission is expected. 
However, GMCs are quite porous to UV radiation (Boisse 1990) 
so we expect to see a rather thick cloud surface, made thicker 
in a low metallicity object like NGC 6822. Cloud cores make 
up only a small fraction of the molecular mass of a galaxy and 
this is particularly clear for NGC 6822 from the weak '^^CO and 
HCN emission (Sect. 19. 21 ). As shown in Sect.|4] we see individual 
GMCs in NGC 6822 (individual because the narrow lines cannot 
come from an accumulation of objects) and the majority of them 



are of order our beam size, i.e. not unreso lved, and similar in 
size to Galactic GMCs. iBendo et all (l2008h show that the PAH 
emission at large (kpc) scales in nearby spirals seems to trace 
the cool diffuse dust responsable for most of the I60 Afm emis- 
sion, thereby tracing the gas mass. [Regan et al.l (l2006l) conclude 
from their observations that the PAH emission at large scales can 
be used to trace the interstellar medium. Thus, from both an ob- 
servational and theoretical point of view, at large scales the PAH 
emission can be used to trace neutral gas. The lower metallicity 
in NGC 6822 is not an issue for our method because we use re- 
gions with little or no star formation and low 70 and 160 micron 
emission, such that little or no molecular gas is expected, in or- 
der to "calibrate" the 8yum emission per H-atom per FUV ratio. 

We compute the emissivity of the PAHs per hydrogen atom 
and per unit of ISRF (traced by the GALEX FUV data) in the 
8jum IRAC band, ((S s^,,,,/ FUV)/N{m))o, in regions far from 
major star forming regions and with low but well-measured H I 
column densities. In these regions we find (A^(H I)/5 8^m)o = 
1.25 X 10^^ H cm 2 MJ y'^ sr, close to the value 1.1 x 10"^^ 
erg s ' sr ' H ' given bv lDraine&Lil(l2007h in their Table 4. 

The value of the interstellar radiation field at each position 
was derived from the far ultr aviolet GALEX ma p and the mean 
value of the ISRF given by ( iDraine et alj|2007l) for the whole 
galaxy. No corrections were made to correct for UV extinction 
by Milky Way dust. The GALEX FUV data were taken from the 
GR5 public release of the MAST archive 

Then, from the 8jum per H-atom emissivity, the 8/zm emis- 
sion, FUV emission and the H I column density at each position, 
the molecular gas column density can be estimated at each posi- 
tion as follows: 

2 ^(H2) = l ^^^^^l,] S ,„n/FUV - N{H I) (6) 

At the 15" resolution used for 8;um and FUV, the atomic 
line data is assumed to be sm ooth thus the origin al 12" x 42" 
HI map is directly subtracted. IDraine &Lil ( I2007h find that the 
emissivity per H-atom and per unit of interstellar radiation field 
is constant over a wide range of values of the ISRF and that 
this should be true irrespective of whether the gas is atomic or 
molecular as long as it is neutral. 

Taking into account the variation of the PAH emissivity with 
the radiation field does not (significantly) change the computed 
masses of individual clouds except in the case of Hub V (cloud 
2 in Table |2]i where the mass is found to be 5 times smaller than 
for the non ISRF corrected case, for an estimated interstellar ra- 
diation field of 50 Habing. 

In extreme radiation fields, PAHs can be destroyed. 
However, in NGC 6822 very little of the dust mass is exposed to 
such fields; according to Draine et al. 2007 Table 5, less than 1% 
of the dust in NGC 6822 is exposed to a high ISRF. Furthermore, 
there is no evidence of PAH destruction through low 8/24 //m 
ratios (cf. Table 2 of Cannon et al. 2006). PAH destruction is 
unlikely to affect our estimate of the H2 mass. 

The spatial correlation between the 8 fim PAH emission and 
the FUV is quite good. The presence of zones with an FUV 
peak but without a PAH peak does not affect our mass estimates 
and the opposite, PAH emission adjacent FUV emission coming 
from far enough away that the division would not affect the same 
pixels, would only cause us to overestimate the H2 mass. 

The corresponding H2 column density map is shown in 
the upper panel of Fig. |7] which shows the H2 distribution in 
NGC 6822 derived as above along with the CO(2-l) integrated 
intensity contours for the part observed by us in CO. 
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This way, we estimate the molecular gas mass within 
our observed zone to be M(H2) ~ 6 x 10^ M0 and about 
M(H2) ~ 1 X 10^ Mq over NGC 6822 as shown in Fig.]?] This 
is 20 - 25 times greate r than using a Galactic value and is in ex- 
cellent agreement with llsraell (1 1 997al) who estimated the total H2 
mass to be about 1.5 x 10^ M0. 

We have also calculated the H2 masses for the individual 
clouds in Table 2 using the S/zm emission and the last col- 
umn gives the corresponding estimated A^(H2)//co(2-i) values. 
The median and average values are 15 - 20 times the Galactic 
value, assuming a line ratio of 0.7 to go from A^(H2)//co(2-i) to 
A^(H2)//co(i-o)- Overall, the agreement between the CO emis- 
sion and the S/zm based vision of where the molecular gas is 
found is quite good, suggesting that CO traces the H2, albeit with 
a much higher A^(H2) / /co(i-O) factor than in the Gal axy. Ga lactic 
CO emission is present towards NGC 6822 (lsraei [T997ah (see 
Sect. [TOb . the H 1 column is of course unaffected but the 8//m 
continuum could be. However, since the morphology of the H2 
column density map resembles NGC 6822 so closely but not the 
local emission, we expect this contribution to be low. 

6.4. Molecular gas mass from 160fim emission 

We can perform the same sort of calculation based on the 160;um 
emission, which we assume to come from dust grains large 
enough that they reach thermal equilibrium and are not af- 
fected by in dividual photons. T he pr operties of the grain s can 
be found in iLi & Draind ( 1200 lb and iDraine & Lil (|2007|) . The 
first step is to measure the dust temperature. Since the work by 
llsrae l' ('1997a'), the Spitzer data for NGC 6822 has become avail- 
able (iCannon et al.in2006i) . extending to longer wavelengths than 
IRAS and thus more sensitive to coo l dust. Using t h e 70;u m and 
160/im Spitzer data, we derive, like ICannon et al.l (l2006l) . dust 
temperatures around 25 K (assuming that the dust cross-section 
varies with A'^). Since cool dust emits very weakly at 70 fim , we 
chose to use the 160;um Spitzer and 100 fim IRAS data to better 
measure the temperature of the cool dust component. However, 
some of the 100/im emission may still come from a warm com- 
ponent, causing an overestimate of the dust temperature of the 
cool component and a correspo nding underestimate of the gas 
mass. NGC 6822 has an SED dCannon etalJl2006l figure 12) 
similar to that of NGC 4414, as measured by the ISO LWSOl 
scan byi rame & Hughes! ([T999h . In their Figure 4, they present 
a breakdown of the dust emission into warm and cool compo- 
nents. The 100 yum emission due to the warm dust is about 8% 
of the total 70 micron emission. In equation|7J we therefore take 
X - 0.08 as our fiducial value but also test;^' - andx = 0.16 to 
measure the effect of an error in;^'. We smoothed the 160/Ym data 
to the resolution of the IRAS 100 fim HiRes maps (1 18"x 102"). 
Assuming a modified grey-body law with a spectral index y6 = 2 
for the dust, we then estimate dust temperatures around 2 1 K with 
Eq.[8]instead of 23K using the 70 and 160 fim data for the same 
regions. Changing from to 0.08 causes a 0.3K change in dust 
temperature. 



Tdust - f 



J 100/im„„; 



(7) 



(8) 



The computed dust temperature map is shown at the right of 
Figure[8] The temperature map does not cover the whole galaxy 
because it was necessary to clip the very low signal-to-noise re- 
gions. Cuts were applied at 4 MJy/sr for both the 100 and 160 



fim maps, leaving the odd shape seen in Figure[8] From the dust 
temperature, the HI column density, and the 160 fim emission, 
assumed optically thin, we can estimate the dust cross-section 
at 160 fim which can then be applied to regions with molecular 
gas. 

Although the map of dust temperature is at IRAS resolution, 
we apply it to the full resolution 160 fim data so that the mor- 
phology is better reproduced. Smoothing does not affect the to- 
tal 160 fim flux. Averaging the HI/160mu ratio map over regions 
with high enough 160 emission that the noise has little effect (4 
MJy/sr) and low enough 160 emission to exclude regions where 
molecular gas will be present (8 MJy/sr), we obtain dust cross- 
section at 160 fim of 0-150 = 5.1 x lO^^^cm^ per H-atom. We 
also varied the threshold values (from 4-8 to 6-12), this lead to 
a variation of sigma of at most 10%. 

ICannon et"al] (l2006h found consistent values of Mdust/MHI 
for the individual regions they observed and some of the scatter 
is certainly attribuable to the molecular gas that they could not 
measure, although in most cases the HI dominates. They found 
a factor 5 lower ratio when they calculated Mdust/MHI for the 
whole galaxy. Using the 70/160 dust temperature like Cannon 
and the 160 micron emission, but only over the area where we 
felt we could reliably estimate the dust temperature (cuts at 1 
MJy/sr at 70 fim and 4 MJy/sr at 160 fim , both smoothed to 
the 160 fim resolution) we find a dust cross-section (equivalent 
to Mdust/MHI r atio) of 4.7 x 10_ff_^cm^ per H-atom with little 
variation, unlike ICannon et aL I (l2006l) . 

Assuming a linear dependence of crigo on metallicity 
(Oxygen abunda nce), a solar oxyg en abundance of 12 + 
log(0/H) ^ 8.67 dAsplund et alj2005h . l2+log(0/H) = 8.11 in 



NGC 6822 dLee et alJl2006l) 7 and a solar metallicity dust cross- 



section of o-i(,o - 2.25 x 10 cm- per H-atom (Table 6 of 
iLi & Draind I2OO lb . we obtain am = 6.3 x 10 -*cm- per H- 
atom for NGC 6822. Thus, our "observational" results are in 
good agreement with model calculations. 

Using this temperature map and the cross-section 0-150 = 
5.1 X 10"-*cm^ per H-atom above, the total H column density 
at each position is 



NHtot = S 



1 



160,„„ oi60//m,T<i 

SO that the H2 column density is simply 

Nh: = ^(Nmot-NHi) 



(9) 



(10) 



The left panel of Fig. [8] shows the total H2 column density de- 
rived in this way using an Nhi map smoothed to the 160yum 40" 
resolution. 

In Table 2 the H2 mass estimates for individual clouds are 
only based on the 8/^m map due to angular resolution - the 160 
fim data is at a resolution larger than the clouds. Tables 3 and 
4, which provide a summary of the molecular gas mass calcula- 
tions, include the 160 fim results. 

It is very difficult to estimate uncertainties for our mass esti- 
mates. Statistical noise related errors, as manifested by the vari- 
ations within "blank" areas of maps, are about 10"'^ in column 
density. Systematic uncertainties are certainly present as well 
so we consider the column density noise in our maps to be 
~ 2 X lO^o H2 cm-2. 
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N(H2) from S^eo with T^ust map (cm ) 
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Fig. 8. H2 column density map estimated from 160/im emission using the dust temperature map shown on the right (left), CO(2-l) 
intensity contours every (O.lSKkm s"') starting from (0.15K km s"'). Temperature map from 160yL(m Spitzer and lOOyumlRAS data 
(right) . The yellow rectangle is the region we have mapped in CO(2-l) with HERA. The beams for the NH2, CO and temperature 
are plotted in white. 



Table 3. Average over the cloud sample of the A^(H2)//co(2-i) 
factor for the different methods used to derive the molecular gas 
mass. 



Table 4. Molecular gas masses, A^(H2)//co conversion factor 
and characteristic time to transform molecular gas into stars de- 
rived from IR emission. 



Method 




(A'(H2)//C0(2-1)) 

cm"2/K km s"' 


cm 


Range 

-^Kkms-' 


total Mh, 
Mo 


map Mh, 
Mo 


N(H2)/lco 

K km s-1 


T 

yr 


Sfjm 


by eye 
CPROPS 


2.1 ±0.8 X 10^' 
3.5 ± 1.4 X 10^' 


0.3 
1.0 


- 4.3 X 10^' 

- 6.4 X 10^' 


Sfim 9.2 X 10^ 
160yum 4.6 X 10^ 


5.7 X 10* 
3.4 X 10« 


5.3 X 10^' 
2.3 X 1021 


2.2 X 10** 

1.3 X 10** 


160/jm with 


by eye 


1.5 ± 0.9 X 10^' 


0.5 


-3.1 X 10^' 


The second column refers to the area mapped in CO. The last column 


Tdust map 


CPROPS 


1.3 ±0.9 X 10^1 


0.4 


- 2.7 X 10^' 


is the inverse of the SFE. 








Mvir 


by eye 


1.6 ±0.8 X 10^' 


0.2 


- 2.8 X 10^' 










CPROPS 


1.1 ±0.8 X 10-' 


0.07 


- 3.5 X 10" 











Averages for the "by eye" and CPROPS methods 
1 1 and 9 clouds. In the case of the virial masses 
in the sample are 12 and 1 1 respectively for the 
methods (see TableO. 



are respectively over 
the number of clouds 
"by eye" and CPROPS 



6.5. The A^(H2)//co(2-i) ratio in NGC 6822 

We have estimates of the CO-to-H2 conversion factor on two 
different scales, at the cloud level and for the whole area mapped 
by HERA in CO(2-l). 

Table [3] presents averages and total ranges of the 
A^(H2)//co(2-i) factor for our sample of clouds for the differ- 
ent methods we have used. The average is over 9 clouds in the 
case of CPROPS and 14 in the case of the "by eye" identifica- 
tion. For the virial mass the average is over our full sample of 1 1 
and 15 clouds using a projected area as defined in Sec. 14.2.11 We 
have included the virial mass as a valid method to estimate cloud 



masses because the size-linewidth relationship we find (see Fig. 
nil is similar to the one for Galactic clouds. This suggests that 
the CO molecules are not greatly photodissociated at the outer 
edge of the clouds and that the cloud size as determined from 
the CO emission is close to the true size of the molecular clouds. 
As might be expected in this case, the masses computed by the 
virial theorem are similar to but usually lower than the ones we 
estimate by the other methods (see Table|2|. 

We find similar values of A^(H2)//co(2-i) for the different 
methods used, around 1.5-2x 10^' cm"^/(K km s"'), 5-10 times 
the » 2 X 10^° cm"^/(K km s~') value for the inn er part of the 



Milky Way in CO( l-O) (e.g. iDickman et al.lll986h . From inter- 
ferometric CO data lBolatto et al. (l2008l) using onlv virial masses 
for a sample of dwarf galaxies of the Local Group find a value of 
the CO-to-H2 factor similar to the Galactic value despite the low 
metallicities. Presumably this is because they detect the dense 
protected parts of bright clouds. 
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There has been evidence for a long time that as the size scale 
increases, at least for low metallicity objects, the A^(H2)//co fac- 
tor increases as well ((kubio et al. 1993). Table|4]shows the mass 
estimates for the area of NGC 6822 mapped by HERA and a 
larger region including the HERA map and which covers vir- 
tually all of the stellar and Ha emission. The A^(H2)//co(2-i) 
factors and characteristic times to transform molecular gas into 
stars are shown in columns 3 and 4 for the different methods 
used to estimate the molecular gas column density. After sub- 
traction of the HI column density, some pixels became nega- 
tive; we have attributed a nil value to the pixels in the NH2 
map where the molecular gas column density had a negative 
value after application of Eq. 4. The A^(H2)//co(2-i) factor for 
the whole of the area mapped in CO(2-l) is slightly larger 
(~ 5 X 10^' cm"^/(K km s"')) than for the individual clouds 
and 20-25 times larger than the standard Galactic value for the 
molecular ring of the Milky Way. This high er valu e is clos e to 
the 4 - 8 xlO^' cm-2/(Kkm s ') found by llsraell (Il997al) for 
scales of ~ 50". We now turn to modeling to better understand 
the relationship between CO emission, H2 mass, and the other 
properties of the clouds. 



7. Modeling the CO emission from NGC 6822 GI\/ICs 

7.1. Description of the models 

We have modeled the structure and the emission of the molec- 
ular clouds in NG C 6822 using CLOUDY version 07.02 
dPerland et al.l Il998l) . CLOUDY computes the spectrum of a 
gaseous nebula using as only inputs the geometry, the gas 
and dust composition and the energy input on the nebula. 
Originally designed to deal with photoionized nebula, CLOUDY 
is now able to perform accurate calculations in the Photo 
Dissociation Region (FDR) and even well into th e molecular 
clouds due to the implementation of the H2 physics ( IShaw et al.l 
I2005h and a network of 1000 reactions involving 68 molecules 
dAbel et al1l2005h . CLOUDY determines steady-state solutions 
for the chemistry and predicts the emission in the fine structure 
lines and CO, which are the main coolants in the FDR and the 
molecular cloud. 

The exact geometry of the regions to be modeled, the po- 
sition of the excitation sources and the shape of the incident 
continuum are not well known. Therefore we have performed 
pl ane-parallel m odels with the incident continuum determined 
bv lBlackl(ll987h for the local interstellar radiation field (ISRF). 
This continuum is not restricted to the FUV range as in many 
classical FDR calculations. To estimate the intensity of the in- 
ci dent continuum we have taken into account that according 
to iDraine et all (l2007l) . the gas in NGC 6822 is exposed to a 
minimum radiation field of order twice the local ISRF. On the 
other hand, the variation in the GALEX FUV observations of 
NGC 6822 is about a factor 50 from the general field to the 
stronger fields around specific positions. Therefore, we have 
computed simulations for scaled versions of the local ISRF 
by factors of 1, 10, and 100. However, most of the clouds in 
our regions of interest are exposed to (GALEX-based) fields of 
less than 10. In addition to the Black (1987) continuum, the 
calculations include the cosmic microwave background and a 
cosmic -ray i onization rate of 2.5 1 0~'^s~' , which is typical of the 
Milky Way (jWilhams et al.|[T998l) but not known for NGC 6822. 
Given that UV photons penetrate further into molecular clouds 
in the low-metallicity environment of NGC 6822 than into their 
Galactic counterparts, the cosmic ray flux is not as critical a pa- 
rameter 
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-n= 10000, G=l - 
■n= 10000. G=l - 



-- 0=10 G=100 

-- G=10 0=100 

0=10 0=100 
I UMIST 



10^ ^ 
Nh2 (cm ^) 

Fig. 9. Model A^(H2)//ci (top) and A^(H2)//co factors as a func- 
tion of visual extinction A,, and H column density into the cloud. 
As the atomic Hydrogen skin is very thin even for the highest 
UV field used here, the clouds are essentially completely molec- 
ular so the depth is shown as H2 depth not H. The jc-axis extends 
beyond standard cloud column densities in order to show how a 
change in cloud structure with respect to Galactic clouds would 
affect the CO intensities. The expected range is roughly indi- 
cated by a box in the lower panel. While the UV field has little 
effect on the Hydrogen, it has a major effect on the C^/C/CO 
transition. 



We have adopted typical gas abundances and dust distribu- 
tions (including FAHs) of the Milky Way interstellar medium 
(see CLOUDY documentation for the exact values). These abun- 
dances are scaled using the metallicity of NGC 6822 (0.3 solar). 
The simulations have been computed for two densities, 1000 
and 10000 mol cm"^, which are representative of the bulk of 
the molecular gas. In addition we have computed simulations 
with two series of rate coefficients - the CLOUDY default and 
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Fig. 10. Line strengths of the major observables C II, O I, C I, and 
CO as a function of visual extinction A,, and H column density 
into the cloud. Lines are as in Fig. 8 -black for density n - 1000, 
red for the nominal « = 10000 with CLOUDY default rates, and 
green for n = 10000 with UMIST rates. The solid, dashed, and 
dotted lines are for respectively UVfieldsofl, 10, and 1 00 times 
the estimated solar neighborhood value. 



the UMIST. The reason for testing two sets of rate coeffic ients 
is that in the comparison of PDR codes (iRollig et al.ll2007|) . the 
C I-CO transition was found to vary significantly (from A,. ~ 2 
to Av ~ 5) depending largely on these coeflicients. This is a crit- 
ical region for low-metallicity molecular clouds because the CO 
emission depends strongly on the fraction of the cloud in which 
CO has formed. 

In summary, we have explored a grid of 12 models (3 fields 
X 2 densities x 2 reaction rates). In all the simulations, a line 
width of 2 km/s is included via the turbulence command. The 
results are summarized in Figs. |9] and [TO] In both figures, the 
simulation results are shown for three radiation fields, solid lines 
for Go = 1, dashed for Go = 10 and dotted for Go = 100 (where 
Go is the intensity of t he far u l traviol et continuum in the cloud 
surface in units of the iHabingl d 19681) field). The red and black 
curves use the CLOUDY default rate coefficients for respectively 
densities of 10"* and 10^ cm"^. The green curve uses the UMIST 
rate coefficients for lO"* cm""*. 

The range of H2 column densities explored goes from ~ 
3 10^' cm"^ to ~ 3 10^^ cm"^. Assuming that the grain abun- 
dance scales with metallicity, those column densities correspond 
to visual extinctions (A,,) from 1 to 10. Therefore, the model sim- 
ulations cover the full range of H2 column densities expected in 
the region of interest, from a few 10^' cm"^ of the extended com- 
ponent derived from the dust emission to ~ 10^^ cm"^ towards 
the center of the molecular clouds detected in CO (see Fig. |7ll. 
The l atter value is a lso typical for molecular clouds in the Milky 
Way ( lLarsonlll981h . although regions of higher column density 
also exist, particularly towards cloud cores and GMC centers. 

Among the goals of the modelling is to see if the A^(H2)//co 
factor is close to our other estimates and to predict testable fluxes 
in other important cooling lines which could then act as further 
diagnostics. Figure|9]shows the N(H2)/Ico and N{H2)/Ici ratios 
as a function of the total column density and visual extinction 
into the cloud for all models. In other words, at a given depth 
into the cloud the plotted N{li2)/Ico and A^(H2)//ci ratios are 
the integrated values from the cloud surface to that cloud depth. 
While, Fig. [To] shows the cumulated intensity for different lines 
as as a function of the total column density into the cloud. 

The column densities shown in Figs. |9] and [TO] are actually 
total (H + H2) but the gas is almost completely molecular due 
to the low radiation fields (the model with a highest atomic H 
fraction has only 3 % of H in the region of interest). 

The C I to CO transition occurs when the lines in Fig. |9] are 
decreasing steeply due to the fact that CO begins to emit strongly 
just after formation. This is clearer in the two lower panels of 
Fig. \W\ which shows how the energy emitted is divided into 
the principal PDR cooHng lines: CO, C I, C II, and O I[63]. The 
models suggests that no CO emission is expected for H2 column 
densities below ~ 2 10^' cm"^. 

7.2. Constraining tlie pliysical conditions and tlie A^(H2)//co 
ratio 

Currently, the spectral data available to study NGC 6822 is lim- 
ited to HI and CO. However, in this section we will try to con- 
strain the physical conditions and the N{H2)/Ico ratio using our 
model calculations and assuming line ratios typical of normal 
galaxies. In doing so, we will also make u se of GALEX mea - 
surement of the far-UV field in NGC 6822 jDraine et al.ll2007h . 
We will also take into account our previous estimations of the 
H2 column density of the CO clouds detected in NGC 6822. 
The exact value is not important since we show below that the 
A^(H2) / Ico conversion factor depends weakly on NH2 in the re- 
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gion of interest. Therefore, the Af(H2)//co conversion factor de- 
rived in this section remains independent of other previous esti- 
mations. 

For normal spiral galaxies, the CII[158] integrated inten- 
sities in energy u nits are typicallv 1000 - 2000 times th e 
CO(l-O) intensity dStacev et al.lll991l: iBraine & Hugheslll999h . 
On the other hand, the OI[63] to CII[1581 ratio is ~ 0.3 
in normal galaxies at large scales (iBraine & HughesI Il999t 
iLord et al.lll996t iMalhotra et al.ll2001h and can be even higher 
when studying smaller scales as the spiral arms of M31 
dRodriguez-Fernandez et al.]l2006h . These are lower limits to the 
O I/C II ratios in the neutral components of the ISM since a frac- 
tion of the observed C II emission may come from ionized gas 
not closely associated with molecular clouds or PDRs. 

The simulations show that the C II flux varies linearly with 
the incident field and depends little on the density for the range 
explored here, as expected. The OI flux depends strongly on 
both density and radiation field and for the low density mod- 
els the flux is well below the large-scale value of 1/3 that of the 
CII typical of normal galaxies. Thus, a density of 10"* cm is 
more appropriate to reproduce the O I/C II ratio. 

Even for large column densities where the CO is formed for 
all the parameters studied, the simulations predict strong varia- 
tions in the CII/CO(1-0) ratio from 100 to 10'* for incident field 
of Go = 1 to 100 respectively. To reproduce ratios of 1000-2000 
typical of normal spiral galaxies a relatively low intensity inci- 
dent field of Go = 10 is favored. Similarly, the model CI/CO(l- 
0) ratio also favors low intensity fields since the predicted ratios 
are high, about 0.5 in K km/s except for G = 1 with a ratio of 
about 0.3 , whereas typical observed values are about 0.2 (e.g. 
iGerin & P hillips 2000) and often less. For lower columns or the 
late CO formation curves, the C I/CO ratio is substantially higher 
(due to the weak CO). 

All together, the comparison of the model prediction with 
the typical line ratios yields a representative density of ~ 10"* 
cm'^* and in c ident field in the range 1-10, in agreement with 
iDraine et al.l (|2007|) and our GALEX estimates for NGC 6822. 
Taking into account those density and incident radiation fields, 
and a range of H2 column densities from 8 10^' cm"^ to 2 10^^ 
cm which comprise the expected H2 column density of the CO 
clouds detected in NGC 6822, a likely value for A^(H2)//co(2-i) 
is in the range ~ 2 - 3 x 10^' cm"^/(K km s"'). These estima- 
tions of the A^(H2)//co(2-i) factor do not depend critically on the 
exact H2 column density of the molecular clouds in NGC 6822 
and agree reasonably well with other estimations. Increasing the 
radiation field further pushes the CO edge of the cloud deeper 
in, probably unrealistically deep unless the molecular clouds in 
NGC 6822 have higher column densities than Galactic GMCs. 
Observing the C II, O I, and C I lines with the Herschel satellite 
would be of great interest because they are sensitive respectively 
to the SFR, the SFR and density, and the cloud depth (the others 
being independent of total cloud column density - cf. Fig.fTOli. 



8. The efficiency of star formation in NGC 6822 



Table 5. Properties of the H II regions Hubble V and Hubble X. 



Using t he Ha luminosity a nd the calibration from iKennicutl 
( Il998l) . ICannon et al.l ( 20061) derive a global star formation rate 



of .0 1 5 M0 vr' ' . Ilsrael et aP (1 1 996h find 0.04 M© yr ' over the 



last 10^ years from bolometric luminosity measurements. Using 
the calibration from iHun ter & Gallagher ( 1986) which supposes 
that close to half of the ionizing photons are lost to dust and a 
Salpeter IMF: 



Name 




<5off'' 


TcO rab 

mK 


Icq 

K km 


Hubble V 


-68.5 


299.0 


500 


0.93 


Hubble X 


55.2 


310.8 


< 85'' 


< 0.19' 



" Offsets with respect to the reference position 
(a„,6„) = (19''44™57^:8,-14°48'11") 
Corresponding to a 3cr level 

Corresponding to a 3<t level and a 3 km s"' fwhm linewidth. 



yields, with L(H„) = 2 x 10^*^ erg s ' dde Blok & Waltej|2006bh 
a star formation rate of 0.026 Mq yr"'; this is the value we will 
use for NGC 6822. 

The star formation efficiency is usually defined as the ratio of the 
star formation rate over the mass of molecular gas available to 
form stars: 



SFE = 



SFR 

mT 



(12) 



M = 5 X IQ-^L{H„)/Lq Mq yi- 



(11) 



The characteristic time to transform molecular gas into star is 
directly r = 1 /SEE. Table |4] shows the values of the characteris- 
tic time T for the different methods used to estimate the molec- 
ular gas column density. The H2 consumption tim es are smaller 
than those found for la rge spirals, about 2 x 10^ yr ( iMurgia et al.l 
l2002l:lKennicuttll998h . such fliatthe SEE in NGC 6822 is higher 
by the same ratio. This is consistent wit h the high star formation 
efficiency found bv lGardan et al.l (l2007h in M 33, a spiral galaxy 
larger than NGC 6822 but an order of magnitude less luminous 
and less massive than the Milky Way. Two other small galax- 
ies, ICIO and NGC 2403 (which is very similar to M33), have 
also b e en found to have high SEEs by respective ly iLerov et all 
(l2006h . iKennicuttI (Il998h and iBUtz et alJ (l2007h but the con- 
straints on the A^(H2)//co conversion are not stringent. 

Going to z ~ 1, the sta r formation rate increases by at 
least one order of magnitude (iMadau et al.l 1 19961 : IWiUdns et al.l 
I2OO8I). Even if the gas fraction was higher in the past, a high 
star formation efficiency has to be introduced in order to ex- 
plain such a wide variation of the star formation rates. Small 
Local Group galaxies such as NGC 6822, M 33, the LMC or 
the SMC share some properties with intermediate redshift ob- 
jects: they are gas rich, have subsolar metallicities and seem to 
exhibit high star formation efficiencies. Low-luminosity NGC 
6822 also shares with early universe objects a high FIR/CO lu- 
minosity ratio. Note that both CO and the dust emission are af- 
fected by metallicity. However, unlike these rare but very lu- 
minous galaxies, NGC 6822 (like the other small local objects 
above) ha s a low SFR, placing it in an em pty region of Figures 
8 and 9 in lSolomon & Vanden Bouj (l2005l) . High spatial resolu- 
tion observations of these local systems may help us understand 
the physics of intermediate redshift galaxies. 

9. Hubble X and Hubble V 

9. 1 . Diffuse CO emission Soutfi of Hubble X 

A major question is whether large quantities of H2 could be 
missed through the use of CO as a tracer of molecular gas. In 
particular, does the presence of the luminous star forming re- 
gion Hubble X (see Table |5] for positions of Hubble V and X) 
with copious Ha, FUV, 160/im and 8//m emission but without a 
CO detection invalidate CO as a tracer of H2 in galaxies like 
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Fig. 11. CO(2-l) and HI spectra of the region south of of Hubble 
X in which CO emission appears to be detected (cf. Fig.|6]l. The 
solid line is the CO(2-l) summed over a 2.57 arcmin^ region 
east of Hubble X, the dotted line is the H I spectrum summed 
over the same region, and the dashed line is a JCMT CO (2-1) 
spectrum taken at RA 19''42"n9 -^.7, Dec -14°50'31". 
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Fig. 12. Spectra of the Hubble V (iHubbld [19251) HII region. 
(dotted line) '^CO divided by 10, (continous line) ^^CO, (black) 
J = 2 -> 1, (red)J= 1 ^ 0. 



NGC 6822? As shown in Fig. |7l for reasonable estimates of 
the radiation field or dust temperature, the dust emission corre- 
sponds to that expected from the atomic component alone. Thus, 
we beUeve that towards the H II region there is in fact very little 
H2, and not just little CO. However, we see evidence for "dif- 
fuse" CO emission from the Hl-rich region near and South of 
Hubble X (see Fig. [nj. 

In Fig.[TT] we have summed the CO and HI spectra over a 
region roughly 50" x 200" in size, yielding an apparent detection 
in CO despite the absence of individual clouds with detectable 
CO emission. Also shown is the result of an observation towards 
19H2'"19 ^7, -14°50'31" made with the JCMT in CO(2-l) at 
21" resolution in 1992, which, although the signal-to-noise ratio 
is low, seems to show CO emission towards roughly the same 
position and at the same velocities. With the resolution available 
in extragalactic observations, "diffuse" signifies that neither spa- 
tially nor spectrally can we distinguish the clouds which, taken 
together over a large area, appear to contribute a detectable CO 
signal. While these clouds could be like Galactic cirrus, our res- 
olution and brightness sensitivity are not sufficient to be sure. 
They are not like either galactic GMCs or the individual clouds 
discussed in Sections H] Using the methods described in Sect. 
6, no H2 is found in this area although the CO luminosity is 
~ lO'^Kkm/spc". 

9.2. The "CO and HCN emission in Hubble V 



Unlike most spiral galaxies (iBraine et al ] Il993h . compar- 
ing the CO(l-O) emission with the CO(2-l) emission con- 
volved to 21" angular resolution yields a line ratio of about 
CO(2 - l)/CO(l - 0) ~ 1.2 in Hub V. This is unlikely to be the 
case over much of NGC 6822 because rather warm and dense gas 
is required, including some optically thin emission. The '^CO 
emission confirms this for Hubble V - the '^CO/'^CO ratio is 15 



in the (1-0) transition but only 5.7 in the (2-1) trans ition, show- 
ing that the higher transition is efficiently excited. Ilsrael et aH 
(2003a) found considerably higher '^CO/'^^CO line ratios with 
the SEST telescope. This is a strong indication for the presence 
of tenuous molecular gas with low CO optiocal depths surround- 
ing the Hubble V cloud, probably similar to the diffuse gas found 
near Hubble X (see Sect. 19711) 

In order to estimate the fraction of dense gas in Hubble V, 
we observed the HCN(l-O) line at 88.6316 GHz. HCN has a 
high dipole moment and thus requires high densities to be ex- 
cited. HCN(l-O) was not detected despite reaching a 3cr noise 
level of 0.035 K km s"', a factor 60 below the CO(l-O) fine. 
Typic al values of CO(1-0) /HCN(1-0) are ~ 50 in gala ctic 
disks dKuno & Nakail 1 1 997L for M 51. lBraine et"al1ll997l for 
NGC 4414. iBrouillet et alf^OOl for M 31), 10 - 20 in nuclei 
dNguven et alJll992t iHenkel et al.lll99lh . and less in ultralumi- 
nous infrared galaxies which are as a result believed to have a 
partic ularly high fraction of their Ht in dense co res jGao et al.l 
I2007h . A gain with the SEST, Ilsrael et aP (l2003ah reports detec- 
tion of the HCO^(l-O) line in Hubble V, at4 times the brightness 
of our hmit to HCN(l-O). While rather extreme, Brouillet et al] 
(I2005h noted an apparent rise in the HCO^/HCN ratio going to- 
wards the outskirts of M31, where the metallicity presumably 
decreases. From the data presented in this paper, despite being 
a major H II region, Hubble V is not particularly rich in dense 
gas. If the CO underestimates the H2 mass in Hubble V, then 
the dense gas fraction in Hubble V is lower than in spiral disks 
because the HCN in dense cores should be less affected by the 
dissociating radiation field than the CO. 



10. Galactic emission 

A factor which greatly complicates studies of NGC 6822 is the 
presence of a Galactic molecular cloud along the line-of-sight 
towards NGC 6822. Some optical/UV emission is absorbed and 
the cloud emits at (at least) FIR wavelengths, making it diffi- 
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Fig. 13. CO(2-l) emission from the local molecular "cirrus" 
cloud along the line of sight to NGC 6822. Indicated on the fig- 
ure are the velocities of the emission at 3 different positions. 



cult to clearly identify the emission coming from NGC 6822. 
This is particularly a question at 160yL(m where the dust tem- 
perature apparently decreases to the West but this is likely due 
to the increasing column density and low temperature of the 
Galactic cloud. In CO, the Galactic cloud can be separated from 
NGC 6822 (see Fig. [TTb so that we know the western part is 
more affected. A velocity gradient from roughly East to West is 
present in the cloud, in addition to the column density gradient. 
For a cloud distance of 100 pc, realistic due to the -20° galactic 
latitude of NGC 6822, the resolution is about 1500 a.u. or 0,007 
pc, making this one of the highest-resolution observations up to 
now of a local cloud. A dedicated study of the large scale struc- 
ture of this Galactic cirrus cloud will be presented by Israel et al. 
(in prep). Figure [13] shows where FIR emission and optical/UV 
absorption are most expected. The decrease in dust temperature 
at the NW corner of the 160yum map in Fig.[8]might be due to the 
cirrus. However, from looking at the region around RA 19:44:32 
Dec -14:45:00, it is clear that any emission from the cirrus falls 
below our intensity cuts, whether at 100 or 160//m. 



11. Conclusions 

From large-scale CO mapping of a very nearby low metallicity 
galaxy, we identify a sample of molecular clouds, most of which 
were also detected by the CPROPS algorithm, allowing an un- 
biased assessment of their properties. The properties of these 
GMCs (size, linewidth, virial mass) appear similar to Galactic 
clouds but the H2 mass per CO luminosity is much higher such 
that we estimate A^(H2)//co(2-i) ~ 2 x 10^' cm"-/(K km s"') 
for the clouds. A variety of methods yield coherent values. Our 
modeling with the CLOUDY program also provides similar ra- 
tios for reasonable parts of parameter space. At large scales 
(~ kpc, well above that of the clouds), we find evidence for a 
higher A^(H2)//co, ~ 4 x 10^' cm~^/(K km s ~') , in agreement 
Wil li work in NGC 6 822 by [Israeli (Il997ah and in the SMC 
by iRubio et"an (Il993h . No evidence for Ht in the absence of 
CO emission was found. The molecular gas masses derived us- 
ing the high N(H2)/Ico values estimated here, coupled with an 
SFR - 0.026 M0 yr"' ±50%, support the idea that molecular gas 
is more quickly cycled into stars in these small low metallicity 
galaxies. In turn, this appears coherent with our image of rapid 
star formation in intermediate redshift galaxies. 
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